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Abstract 

We discuss the transport of matter waves in low-dimensional waveg- 
uides. Due to scattering from uncontrollable noise fields, the spatial co- 
herence gets reduced and eventually lost. We develop a description of this 
decoherence process in terms of transport equations for the atomic Wigner 
function. We outline its derivation and discuss the special case of white 
noise where an analytical solution can be found. 



Introduction 

We discuss in this contribution the transport of atomic matter waves in a low- 
dimensional waveguide. Such structures may be created close to solid substrates 
using electro-magnetic fields: the magnetic field of a current-carrying wire com- 
bined with a homogeneous bias field, e.g., gives rise to a linear waveguide [HI 0, H]. 
Planar waveguides may be constructed with repulsive magnetic [§[] or optical 
fields that 'coat' the substrate surface. The atomic motion is characterised by 
bound vibrations in the 'transverse' direction(s) and an essentially free motion 
in the 'longitudinal' direction(s) along the waveguide axis (plane), respectively. 
Although direct contact with the substrate is avoided by the shielding potential, 
the atoms feel its presence through enhanced electromagnetic field fluctuations 
that 'leak' out of the thermal solid, typically held at room temperature. We have 
shown elsewhere that these thermal near fields are characterised by a fluctua- 
tion spectrum exceeding by orders of magnitude the usual blackbody radiation 
[H H) H 3 • The scattering of the atoms off the near field fluctuations occurs at a 
rate that may be calculated using Fermi's Golden Rule. The consequences of mul- 
tiple scattering is conveniently described by a transport equation that combines 
in a self-consistent way both ballistic motion and scattering. 
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The purpose of this contribution is to outline a derivation of this transport 
equation. The status of this equation is similar to that of the quantum-optical 
master equations allowing to describe the evolution of the reduced density matrix 
of an atomic system, on a time scale large compared to the correlation time of the 
reservoir the system is coupled to, typically the vacuum radiation field. In the 
case of transport in waveguides, we face both temporal and spatial dynamics and 
therefore restrict our attention to scales large compared to the correlation time 
and length of a fluctuating noise potential. Our analysis uses a multiple scale 
expansion adapted from fllQf . Similar to the quantum-optical case, we make an 
expansion in the perturbing potential to second order. In the resulting transport 
equation, the noise is thus characterised by its second-order correlation functions 
or, equivalently, its spectral density. In the case of white noise, the transport 
equation can be explicitly solved. We have shown elsewhere || that this approxi- 
mation holds quite well for thermal near field fluctuations. For technical noise, it 
also holds when the noise spectrum is flat on a frequency scale roughly set by the 
'longitudinal' temperature of the atoms in the waveguide. The explicit solution 
yields an estimate for the spatial coherence of the guided matter waves as a func- 
tion of time. The paper concludes with some remarks on the limits of validity of 
the present transport theory. It cannot describe, e.g., Anderson localisation in 
one dimension JTX|] because on the coarser spatial scale of the transport equation, 
the scattering from the noise field is assumed to take place locally; interferences 
between different scattering sequences are not taken into account. Decoherence 
in 'curved' or 'split' waveguides also needs a refined theory because of the cross- 
coupling between the transverse and longitudinal degrees of freedom, the former 
being 'frozen out' in our framework. 



1 Statistical matter wave optics 

The simplest model for atom transport in a low-dimensional waveguide is based 
on the Schrodinger equation 

h 2 

i^(x, t) = -^V 2 ^ + ^(x, t)il> (1) 

The coordinate x describes the motion in the free waveguide directions. The 
transverse motion is 'frozen out' by assuming that the atom is cooled to the 
transverse ground state. Atom-atom interactions are neglected, too. V(x., t) is 
the noise potential: for a magnetic waveguide, e.g., it is given by 

V{x,t) = (s\p-B(x,t)\s), (2) 

where \s) is the trapped internal state of the atom (we neglect spin-changing 
processes), and B(x, t) is the thermal magnetic field. The noise potential is a 
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statistical quantity with zero mean and second-order correlation function 

C v (s,t) = (V(x + s,t + r)V(x,t)), (3) 

where the average is taken over the realisations of the noise potential. We assume 
a statistically homogeneous noise, the correlation function being independent of x 
and t. As a function of the separation s, thermal magnetic fields are correlated on 
a length scale l c given approximately by the distance d between the waveguide axis 
and the solid substrate f|. This estimate is valid as long as the wavelength 2nc/u 
corresponding to the noise frequency u is large compared to d: for micrometre- 
sized waveguide structures, this means frequencies below the optical range. The 
relevant frequencies of the noise will be identified below and turn out to be much 
smaller than this. 

The coherence properties of the guided matter waves are characterised by the 
noise-averaged coherence function (the time dependence is suppressed for clarity) 

p(x;s) = <^*(x-fs)^(x + §s)>. (4) 

In complete analogy to quantum-optical master equations, this coherence function 
may be regarded as the reduced density matrix of the atomic ensemble, when the 
degrees of freedom of the noise are traced over. The Wigner function gives a 
convenient representation of the coherence function: 

where D is the waveguide dimension. This representation allows to make a link 
to classical kinetic theory: W(x, p) may be viewed as a quasi-probability in phase 
space. For example, the spatial density n(x) and the current density j(x) of the 
atoms are given by 

n(x) = Jd D pW(x,p) (6) 

j(x) = fd D p±W{x,j>) (7) 
J m 

We also obtain information about the spatial coherence: the spatially averaged 
coherence function T(s, t), for example, is related to the Wigner function by 

r(s,t) = y d D £ p(x; s, t) (8) 
= J d D xd D pe [p - s/h W{x,p,t) (9) 

In the next section, we outline a derivation of a closed equation for the Wigner 
function in terms of the noise correlation function. 
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2 Transport equation 



Details of the derivation of the transport equation may be found in the ap- 
pendix |A[ We quote here only the main assumptions underlying the theory. 

(i) The noise potential is supposed to be weak so that a perturbative analysis 
is possible. As in quantum-optical master equations, a closed equation is 
found when the expansion is pushed to second order in the perturbation. 

(ii) The scale l c over which the noise is spatially correlated is assumed to be 
small compared to the characteristic scale of variation of the Wigner func- 
tion. This implies a separation of the dynamics on short and large spatial 
scales, the dynamics on the large scale being 'enslaved' by certain averages 
over the short scale. Similarly, we assume that the potential fluctuates 
rapidly on the time scale for the evolution of the Wigner function. These 
assumptions correspond to the Markov approximation of quantum optics, 
where the master equation is valid on a coarse-grained time scale. 

The derivation of the master equation is based on a multiple scale expansion. 
Functions /(x) of the spatial coordinate are thus written in the form 

/(x) = /(X,£) (10) 

where X gives the 'slow' variation and the dimensionless variable £ = x/i c gives 
the 'rapid' variation on the scale of the noise correlation length l c . Spatial gradi- 
ents are thus expanded using 

V X = V X + ^V^ (11) 

By construction, the first term is much smaller than the second one. Finally, the 
Wigner function is expanded as 

W(x, p, t) = W (X, p, t) + tj^W^X, & p, t) + 0( V ) (12) 

where 77 <C 1 is the ratio between the correlation length l c and a 'macroscopic' 
scale on which the coordinate X varies. The expansion allows to prove self- 
consistently that the zeroth order approximation W does not depend on the 
short scale and to fix the exponent 1/2 for the first order correction. 

The resulting transport equation specifies the evolution of the Wigner function 
Wo. Dropping the subscript 0, it reads 

(3c + £-V x W(x,p)= (13) 
J dV S v (p' - p, E p , - E p ) [W (x, p') - W(x, p)] , 
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where Sy, the spectral density of the noise, is essentially the spatial and time 
Fourier transform of the noise correlation function 

s >-(^ A£ )^/SF c " (s - r ) e ^ a&,/s - (14 > 

The left hand side of the transport equation gives the free ballistic motion of the 
atoms in the waveguide. If an external force were applied, an additional term 
F • V p would appear. The right hand side describes the scattering from the noise 
potential. E p = p 2 /2m is the de Broglie dispersion relation for matter waves. 
We observe that scattering processes p — > p' occur at a rate given by the noise 
spectrum at the Bohr frequency (E p — E p >)/h. If the potential noise is static (as 
would be the case for a 'rough potential'), then its spectral density is proportional 
to S(AE), and energy is conserved. If we are interested in the scattering between 
guided momentum states, then the initial and final energies E p , E p i are typically 
of the order of the (longitudinal) temperature kT of the ensemble. The relevant 
frequencies in the noise spectral density are thus comparable to kT/%. 



3 Results 
3.1 White noise 

White noise is characterised by a constant spectral density, i.e., the noise spec- 
trum SV(q, AE) is independent of AE. Equivalently, the noise correlation is 
^-correlated in time: 

C v (s,t)=B v (s)5(t). (15) 
The integration over the momentum p' in (O) is now not restricted by energy 



conservation, and the right hand side of the transport equation becomes a con- 
volution. One therefore obtains a simple solution using Fourier transforms. De- 
noting k (dimension: wavevector) and s (dim.: length) the Fourier variables 
conjugate to x and p, we find the equation 

(d t + — - V s ) W(k, s) = - 7 (sMk, s). (16) 
v m J 

where we have introduced the rate 

7 (s) = ^ (*v(0) - B v (s)) . (17) 

Eq.(fl6|) is easily solved using the method of characteristics, using s — hkt/m as 
a new variable. One finds 

W(k,s;t) = Wi(k,s- hkt/m) x 

x exp [- f dt'j(s - hkt'/m)] , (18) 
Jo 
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where W^(k, s) is the Wigner function at t = 0. 

We observe in particular that the spatially averaged coherence function (§) 
shows an exponential decay as time increases: 

r(s;t) = r j (s)exp[- 7 (s)t]. (19) 

We can thus give a physical meaning to the quantity 7(s): it is the rate at which 
two points in the matter wave field, that are separated by a distance s, lose their 
mutual coherence. This rate saturates to 7 = 7(0x3) = B v (0)/h 2 for distances s ^> 
l c large compared to the correlation length of the noise field (the correlation By(s) 
then vanishes). This saturation has been discussed, e.g., in fl2]] . As shown in 0, 



the rate 7 is equal to the total scattering rate from the noise potential, as obtained 
from Fermi's Golden Rule. For distances smaller than l c , the decoherence rate 
7 (s) decreases since the two points of the matter wave field 'see' essentially the 
same noise potential. The exact solution ( |T9"D thus implies that after a time of 
the order of the scattering time I/7, the spatial coherence of the atomic ensemble 
has been reduced to the correlation length l c . The estimates given in |J imply a 
time scale of the order of a fraction of a second for waveguides at a micrometre 
distance from a (bulk) metallic substrate. Significant improvements can be made 
using thin metallic layers or wires, nonconducting materials or by mounting the 
waveguide at a larger distance from the substrate |§ . 

At timescales longer than the scattering time I/7, the spatial coherence length 
of the atoms decreases more slowly, approximately as l c /y/yt @. This is due to 
a diffusive increase of the width of the atomic momentum distribution, with a 
diffusion constant of the order of D = h 2r y/l 2 . This constant is in agreement with 
a random walk in momentum space: for each scattering time I/7, the atoms 
absorb a momentum q c = H/l c from the noise potential. The momentum step q c 
follows from the fact that the noise potential is smooth on scales smaller than l c , 
its Fourier transform therefore contains momenta up to h/l c . 



3.2 Fokker- Planck equation 

The momentum diffusion estimate given above can also be retrieved from the 
transport equation, making an expansion of the Wigner distribution as a function 
of momentum. We assume that the typical momentum transfer q c absorbed from 
the noise is small compared to the scale of variation of the Wigner distribution, 
and expand the latter to second order. This manipulation casts the transport 
equation into a Fokker-Planck form 

(d t + ^ • V x + F dr (p) • V p ) W(x, P ) = (20) 
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where the drift force and the diffusion coefficient are given by 



Fdr(p) = - 

A;(P) = / 



/ 



qqS v (q,E ] 



p+q 




(21) 



d qqiqjSvfaEp+q-Ep). 



(22) 



In the special case of white noise, the p-dependence of these quantities drops out. 
Also the drift force is then zero because the noise correlation function is real and 
the spectrum SV(q) even in q. Since q c gives the width of the spectrum, the 
diffusion coefficient turns out to be of order q 2, y, as estimated before. 

Casting the transport equation into Fokker-Planck form, one can easily take 
into account the scattering from the noise field in (classical) Monte Carlo sim- 
ulations of the atomic motion: one simply has to add a random force whose 
correlation is given by the diffusion coefficient. 

We note, however, that the Fokker-Planck equation cannot capture the initial 
stage of the decoherence process, starting from a wave field that is coherent 
over distances larger than the correlation length l c . Indeed, it may be shown 
(neglecting the p • V x term and the drift force, assuming an isotropic diffusion 
tensor for simplicity) that (|20| ) yields a spatially averaged coherence function 



This result implies a decoherence rate proportional to s 2 without saturation. It 
is hence valid only at large times (compared to the scattering time I/7) where 
the exponentials in both solutions (|l^, ^3|) are essentially zero for s > l c . 

4 Concluding remarks 

We have given an outline of a transport theory for dilute atomic gases trapped 
in low-dimensional waveguides. This theory allows to follow the evolution of 
the atomic phase-space distribution (more precisely, the atomic, noise-averaged 
Wigner function) when the atoms are subject to a noise potential with fluctua- 
tions in space and time. The spatial coherence of the gas can be tracked over 
temporal and spatial scales larger than the correlation scale of the noise, in a man- 
ner similar to the master equations of quantum optics. We have given explicit 
results in the case of white noise, highlighting spatial decoherence and momentum 
diffusion. 

The transport equation has to be taken with care for strong noise poten- 
tials because its derivation is based on second-order perturbation theory. It is 
certainly not valid when the 'mean free path' ~ v/j (v is a typical velocity of 
the gas) is smaller than the noise correlation length l c because then the Wigner 
distribution changes significantly over a small spatial scale. (In technical terms, 
the approximation of a local scattering kernel in ( |T3"D is no longer appropriate.) 




(23) 
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Also, the theory cannot describe Anderson localisation in ID waveguides with 
static noise [|ll|]. This can be seen by working out the scattering kernel with 
S v (q,AE) = S v (q)6(AEy. 

2m J dp' S v (p' - p) 5(p' 2 - p 2 ) [W(x,p') - W(x,p)] 

= mSv ^ [W(x, -p) - W(x, p)] . (24) 
p 

We find a divergence of the scattering rate at p — > since the spectrum Sy{2p) 
is finite in this limit. The one-dimensional, static case therefore merits further 
investigation. We also mention that is has been found recently that Anderson 
localisation is destroyed when time-dependent fluctuations are superimposed on 



the static disorder [13, |l4j]. In this context, transport (or master) equations 



similar to our approach have been used. 

Acknowledgements. We thank S. A. Gardiner, S. Potting, M. Wilkens, and 
P. Zoller for constructive discussions. Continuous support from M. Wilkens is 
gratefully acknowledged. 

A Multiple scale derivation of the transport equa- 
tion 

The Schrodinger equation ([!]) gives the following equation for the Wigner function 
(<9 t + p-V x )W(x,p) = (25) 

4 / (l^ (q ' t] eiqx K (x ' P + ^ q) _ ^ (X ' P " ^ q) . 

where V(q,t) is the spatial Fourier transform of the noise potential. Since this 
potential is assumed weak and varies on a scale given by the correlation length 
l c , we introduce the following scaling 

V{q c u, t)= J d D x e- i(?cU X /V(x, t) = lfv P V(u, t) (26) 

where q c = h/l c is the typical momentum width of V"(q, t) and u is a dimensionless 
vector. The parameter rj is given by the ratio between the small scale l c and the 
'macroscopic' scale of the position distribution, the (positive) exponent f3 remains 
to be determined. We assume 77 <d 1 and make the multiple scale expansion ([12]) 
for the Wigner function. Using the expansion ([□]) for the spatial gradient, we 
get 



(W + rj a W 1 )= (27) 
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~ J ^d^(u, t) e iu ^ [W(x, p + g c u/2) - W(x, P - ? c u/2)] 

We now take the limit 77 — > 0, / c — > at fixed g c . The most divergent term 
on the left hand side is the one with (l// c )V^W / o- It could only be balanced 
with a term on the right hand side involving Wq, but due to the small factor 
t/ 3 , this term cannot have the same order of magnitude. We must therefore 
require that (l// c )V^W / vanishes individually: the zeroth order Wigner function 
is independent of the short scale variable £. 

The next terms on the left hand side contain (77° /Z c )V^Wi and Vxl^o, while 
on the right hand side the leading order is (t]^/I c )Wq. We look for a connection 
between Wq and Wi, and therefore, the left hand W\ term must be more divergent 
than the W term. This is the case if 7] a O(l/l c ) > 0(l/X) ~ r)0(l/l c ). We thus 
conclude that a < 1. Comparing powers of r\ on the left and right hand side, 
we find a = (3, since the vector u and the scaled distance £ are of order unity. 
Therefore we get the equation 

(l c d t + --^)w 1 (X,tp)= (28) 
V m J 

i f (i^?y * 

/ 7^w V ^ *) ^ t^X, P + & u / 2 ) " ^(X, P " & u /2)] 

In the exponential, only the short length scale £ = x// c occurs. We thus find 
that the large scale variable X is a parameter in this equation, and get a solution 
via Fourier transforms with respect to £ and t. In the spirit of the Markov 
approximation, we take the slowly varying Wo (as a function of time) out of the 
time integral 

/°° 
dt e lu)t V(u, t) W (. ..,*)« W (. . . , t)V[u, uo] (29) 
-00 

where V [u, u] denotes the double space and time Fourier transform of the poten- 
tial. We note n, u the conjugate variables for the spatial Fourier transform and 
find the following solution for the first order Wigner function 

W 1 (X,£,p)= (30) 

— hr t^tj • 7 — V^n" (Wo X > P + & K / 2 ~ W ° X ' P " & K / 2 

g c 7 27r 7 (27r) w ik • p/m — u c u; + 

The +0 prescription in the denominator is related to causality: it ensures that 
the poles in the complex oplane are moved into the lower half plane, avoiding a 
blow-up of W 1 . 

This result will be inserted into the next order equation that also links W to 

U + p . v x ) w = 

V m J 
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\r\ 2a r d D u 



/CI 11 >> 
— ^V(u, t) e iu « [W X (K, p + g c u/2) - W^X, £, u - g c u/2)] 



g c / c 7 (2tt 

Note that this equation is scaled consistently if 0(1/ X) ~ r/ 2o C(l// c ) = if^OQ/X). 
This determines the exponent a = |. The result is an equation for Wo only. We 
take the statistical average and make the factorisation 

(V(u, t) V[k, lo) W (X, p)) = (V(u, f) a;]) W (X, p). (31) 

This may be justified heuristically as follows: it seems reasonable that the statis- 
tical average can also be performed via 'spatial coarse graining', i.e., taking an 
average over the small-scale fluctuations of the medium. This is precisely the pic- 
ture behind transport theory: the individual scattering events are not resolved 
but only the behaviour of the matter wave on larger scales. The lowest order 
Wigner function Wo may be taken out of the coarse grain average because it does 
not depend on the short scale £ by construction. 

Finally, we introduce the spectral density S(u,u) of the (scaled) noise poten- 
tial 

(V(u, t) V[k, u]) = (2tt) d S(u, u) e iuJt 5(u + k) (32) 

This allows to perform the integration over k when (|30|) is inserted into fl31|) . 
The result still contains a frequency integral where denominators of the following 
form appear 

1 ~ iqc (33) 



i(u/m) • (p + q c u/2) - il c u + E p+qcll - E p - Tuo - iO 

A second term contains the sign-reversed energy difference. These denominators 
ensure that the kinetic energy change occurring in the scattering is compensated 
by a 'quantum' Tko from the noise potential. 

We write the denominators (|33"D as a 5-function plus a principal part. For the 
classical noise potential considered here, the power spectrum S(u, uj) is even in 
u, so that the (^-functions combine and the principal parts drop out. We finally 
get 

d t + ~- V x ) W = (34) 
m J 

¥ J (27p 5(u ' AE/h) [iy ° (x ' p + qcU) ~ Wo(x ' p)] 

where AE = E p+QcU — E p . It is easily checked that this is the transport equa- 
tion (|T3|), taking into account the relation between the scaled and non-scaled 
noise spectra 

^S v (u, AE/h) = S v (q c u, AE) (35) 
n 



that follows from (13) and (|2B|). 
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